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A  NOT*^.  ON  THE  /APPLICATION  OP  SCH''INr7",R»S  VARIATIONAL 
PRINCIPL"  TO  r^IR/^C'S  EQUATION  OF  T^^":  TL'ICTRON. 
By  H,  E.  I.ioses 

1.   The  Variational  Principle. 

Schwinger's  variational  principle  has  been  used  for 
a  wide  variety  of  problems  involving  wave  motion  in  which 
it  is  desired  to  find  the  amplitude  of  a  scattered  wave 
in  term.s  of  the  incident  wave,"   Schwinger's  method  makes 
use  of  the  fact  that  the  amplitude  of  the  scattered  wave 
satisfies  a  variational  principle.   We  shall  indicate 
this  variational  principle  briefly. 

Let  us  consider  a  vector  space.   ''Ve  may  define  two 
different  inner  products,  the  Hermitian  and  symmetric 
inner  product,  in  this  vector  space.   The  Hermitian  inner 
product  (a,b)  of  two  vectors  a  and  b  is  defined  by  the 
condition  that 

(1)  (a,b)  =  (b,a) 

where  the  bar  indicates  the  complex  conjugate.   The  sym- 
metric inner  product  is  defined  by  the  condition 

(2)  (8,b)  =  (b,a)   . 


See,  for  example,  Schwinger's  unpublished  nuclear  phy- 
sics notes,  or  the  lectures  of  N.  Marcuvit?  in  the  notes 
"Recent  Developments  in  the  Theory  of  'Vave  Propagation", 
Inst.  Aop.  Math,  and  Mech, ,  N.Y.U.,  194S-50.  The  point  of 
view  of  the  present  note  is  cl©se  to  that  of  the  latter 
set  of  notes. 


2, 


If  we  w«rk  in  a  vector  space  with  a  Hermitian  inner  prod- 
uct, an  operator  K  is  said  to  be  Hermitian  if 

(3)  (a, Kb)  =  (Ka,b)   . 

On  the  other  hand,  if  we  consider  a  space  which  has  a 
symmetric  inner  product,  an  operator  K  is  said  to  be  sym- 
metric if 

(4)  (a,Kb)  =  (Ka,b)   . 

It  must  be  remem.bered  that  thoufrh  equations  (3)  and  (4) 
hava  the  same  appearsnc^ ,  they  refer  to  different  inner 
products. 

Let  us  consider  a  space  with  a  Hermitian  inner  prod- 
uct.  Let  us  also  consider  the  equation 

(5)  a  =  Ky 

for  the  unknown  vector  y  where  K  is  a  Hermitian  operator 
and  a  is  a  given  vector,   ''e  define  the  real  niimber  A.  by 

(The  second  equation  follows  from  equation  (5)  tof?ether 
with  the  fact  that  K  is  Hermitian.)  I^lso,   we  define  the 
functional   ^\v[  by 

(7)  XW?   =-^Jj4Kv) 

i  )  (v,8)(a,v) 

It  is  easily  seen  that 

(3)  \\7]      =  X   . 

The  variational  principle  used  by  Schwinfrer  states  that 

X  ^v\  has  PS  its stationary  value  the  value  X.   Thus  to 

find  X  one  need  only  find  the  stationary  value  of  X^vj. 
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Prom  (8)  it  is  peen  that  if  v  ppproxinates  y  to  the 
first  order  then  ^W?  approximates  ^\j\    =  ^   to  the 
second  order.   This  fact  may  be  used  to  find  X  to  a  high 
decree  of  accuracy. 

A  similar  variational  principle  is  valid  if  we  work 
in  a  space  ™ith  a  symmetric  inner  product,  and  if  the 
operator  K  is  symm.etric.   In  this  case,  however,  X  Is 
generally  complex. 

For  the  purpose  of  treatment  of  scattering  problems, 
a  generalization  of  the  above  variational  principle  is 
used.   Let  us  again  consider  a  vector  space  with  a  Her- 
mitian  inner  product  and  consider  s    pair  of  equations 

(  a  =  Ky 

(9)  ) 

/  a'  =  K'y' 

where  K  and  K'  are  Hermitian  ad.ioint  operators  which  by 
definition  satisfy  the  condition 

(10)  (K'u,v)  =  (u,Kv) 

for  any  t'^'O  vectors  u  and  v.   If  K'  =  K,  K  is  a  Hermitian 
operator.   It  can  be  sho"'n  from.  (10)  that 

(11)  (a',y)  =  (y',a) 

which  is  called  the  reciprocity  thoorem.   Let  us  define 
a  num.ber  X  by 

and  ths  functional  X\v,v'5  by 

MT^      'V(    t?       (v'  .Kv)        (PvW  .v) 
^"^       n^'^'5  =   (v',aj(a',v)  -=    (v',a)(a'',v) 
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po  that 

(14)  '^\j,Y'l    ="^   . 

It  CPU   be  shown  that  ^7V,v't  is  stationary  for  independent 
variations  of  v  and  v'  about  the  values  y  and  y'  respec- 
tively and  that,  therefore,  from  (14),  the  stationary 
value  of  ^{v,v'r  is  X. 

A  similar  statement  holds  if  the  symmetric  inner 
product  is  used.   In  this  case  K'  is  said  to  be  the  sym- 
metric adioint  operator  of  K  if  K'  and  K  are  related  by 
(10). 

2 ,      Schwiirrp;er '  s  Application  of  the  Variational  Principle 
to  Scattering  Problems* 
It  will  be  convenient  to  introduce  a  notation  for 
operators.   The  elem.ents  of  the  vector  space  in  which  we 
are  interested  will  be  functions  of  the  coordinates,  col- 
lectively denoted  by  the  vector  x,  and  other  variables, 
collectively  denoted  by  y»   '^ri  operator  A  "fhich  operates 

X  V 

on  functions  of  these  variables  will  be  denoted  by  A  »', 
Hence  if  f(x,Y)  is  an  element  of  the  vector  space,  the 
vector  ivhich  results  when  the  operator  A  operates  on  it 
is  A  »'^f(x,Y).   Frequently  it  is  useful  to  consider  oper- 
ators as  being  integral  operators  with  suitable  kernels 
which  may  be  symbolic  functions,  such  as  the  6-functions 
of  Dirac.    e  shall  define  the  kernel  A(x,y>x',y')  ^7 

(15)  A^*"^  f(x,Y)  =  y  A(x,y;x',y'  )  f(x',Y')dx'  dY' 


/  - 
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the  integration  being  taken  over  the  entire  ranp-e  of 
values  of  x  and  y.   The  integration  over  y  is  to  be  in- 
terpreted as  a  sunnation,  if  y  takes  on  discreet  values 
only.   The  condition  that  the  operator  A  is  symmetric 
can  be  expressed  by  means  of  the  kernel  as 

(16)  A(x,Yjx',Y')  =  A(x',y'',^,y)   • 

Likewise  the  condition  that  A  is  a  Hermitian  operator  can 
be  expressed  by  the  relation 

(17)  A(x,yVx',y')  =  A(x' ,y' jx,y)   . 

In  scattering  problems  one  is  interested  in  solving 
differential  equations  of  the  form 

(18)  h'''^X(x,y-,E)  =eX(x,y;E) 

subject  to  boundary  conditions  to  be  discussed  shortly, 

X  V 

Here  H  ''  is  a  Hermitisn  differential  operator,  E  is  a 
real  number  ^vhich  is  considered  fixed,  and  %.{y.,Y'y^)    is 
a  function  of  the  coordinates  denoted  by  the  vector  x  and 
other  variables  collectively  denoted  by  y«   Oi^^  may  re- 
gard the  problem  of  finding  the  solution  %  (x,YiE)  of 

(18)  as  that  of  findinp:  the  eigenfunctions  for  a  fixed 
eigenvalue  E, 

In  scattering  problems  the  operator  H^**''  is  given 
as 

(19)  H^'Y  =  H^'"^  +  q^'"^ 

X  Y 

where  H  *'  is  a  function  of  derivatives  with  respect  to 

the  coordinates  and  can  be  written  E^*"^   =  h''^(t^)  .  and 

O      ©  ax   ' 
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q  '*^  is  a  function  of  x,  which  we  may  write  as 
q''^*'^  =  q'^(x),  which  will  be  assumed  to  vanish  for 
I xl  >  r  . 

This  last  condition  on  q'^Cx)  is  somewhat  severe  ??nd 
for  many  purposes  may  be  replaced  by  the  less  severe  one 
that  q^(x)  approach  rero  sufficiently  rapidly  for  large 
values  of  | x | . 

In  this  manner  we  have  placed  some  restrictions  on 
the  nature  of  the  operator  H  »'  insofar  as  it  operates 
with  respect  to  the  variable  x. 

It  is  now  possible  to  specify  boundary  conditions 
on  the  functions  T'J^^YjE)"   '^he  boundary  conditions  are 
snecified  in  the  following  way:   the  solution  "/^(XjYJ^  ) 
of  the  differential  equation  (18)  is  to  be  such  that  it 
can  be  written  in  the  form 

(20)     'XU,y;E)  =Xin(x,YiE)  +7.  sc(x,Y-,E) 

where  the  function  'Y^  .    (x,y;E)  ("'hich  is  called  the 
incident  wave)  is  a  p-iven  solution  of  the  differential 
equation 

(SI)    Hl'^-/^^^{x,r,^)   =  E  X  i^(-^,Y;E) 

and  where  the  function  7(    (x.YiE)  (which  is  called  the 

'  ^  3  C 

scattered  wave)  is  to  be  found  so  as  to  behave  like  an 
out^Toing  spherical  wave  for  large  values  of  |x|.   These 
boundary  conditions  are  sufficient  to  determine  "X  (x,y>E) 
uniquely.   For  many  problems  in  which  E  is  a  multiple 
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eiffenvalxie  of  the  operator  H^*"^  it  is  useful  to  label  the 
"  .  ■        o 

functions  /.  (x,y'>S)  by  an  additional  set  of  variables. 
The  function  /   (x,Y>E)  satisfies  the  equation 


(22)  (E  -  H^'Y)Y   (x^Y-,E)  =  q'''Y7v(x,Y;S)   . 

Using  an  appropriate  inverse  operator   [E  -  H  »'^]"  ,  equa- 
tion (22)  is  seen  to  be  equivalent  to  the  equation 

(23)  7v3c(x,r-,E)  =  [E  -  hJ*Y]-^  q'''^)(ix,Y;E)   . 

The  inverse  operator   [E  -  H  »']     is  to  be  expressed  in 

o 

terms  of  an  integral  operator  using  for  the  kernel  an 
influence  function  which  gives   /   '^^'^^  proper  behavior 
for  large  values  of  |x|. 

Having  chosen  such  an  Influence  function  wo  desire 
to  express  Y^       for  large  values  of  |x|  as  a  sum  of  terms 
of  the  form  9. (x,y^E)T. (E )  ,  i.e.. 


(24)     7-3c(x,Y-,E)  =  Z  0^(^,Y;E)  T^(E) 

I X  I  ->  GO     , 

where  Q.(x,y>E)  are  appropriately  chosen  spherical  waves 
and  the  functions  T.(E)  are  "amplitudes"  which  are  inner 
products  to  which  one  can  apply  the  second  variational 
principle  discussed  in  Part  I  of  the  present  note. 

The  function  y^(x,YjE)  is  identified  with  the  vector 
y  of  equation  (9),  the  operator  K  is  identified  i.vith  the 
operator   q  '"^  -  q'^*'^[E  -  H^»'^]  q^*"^  ,  and  the  vector  a 
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is  identified  "ith  the  function  q^*"'''/  .„(x,YiE).   The 

in   '  ' 

first  of  equations  (9)  is  merely  a  form  of  equation  (23). 
Tlie  operator  K'  is  tpken  as  the  symmetric  adjoint  or  Her- 
mltinn  adjoint  operator  of  K,   The  vector  a'  will  be  so 
chosen  that  the  fiinction  T.  (E)  is  the  inner  product 
(a',y).   Thus  [T^(E)]~-^  is  identified  with  X.   The  vector 
y'  is  then  tslren  as  the  solution  of  the  second  of  equa- 
tions (9),   It  is  seen  that  the  manner  of  decomposition 
of  /^        into  spherical  wpves  R-iven  by  (24)  determines 

3  C 

how  a'  shall  be  chosen. 


3 .   Application  of  the  Variational  Principle  to  the 
ScatterinfT  of  the  Dirac  Electron. 
A .   Dirac 's  "Equation  for  the  Electron . 

The  above  concepts  can  be  illustrated  by  considering 
the  scattorinrt  of  the  Dirac  electron.   Schwinger  has 
treated  the  problem  Bonewhat  in  his  nuclear  physics  notes. 
However,  he  doec  not  use  an  explicit  form  for  the  inverse 
operator  [E  -  rT*^]"    ,   In  the  present  note  we  shall  make 
use  of  such  an  explicit  form  in  terms  of  an  integral 
operator  and  show  v/hat  quantities  are  to  be  identified 
with  a'  pnd  y* , 

In  Dirac' s  theory  of  the  electron,  the  elements  of 
the  vector  space  are  functions   f(x,Y)  of  the  coordin- 
ates denoted  collectively  by  the  vector  x,  p.nd    of  a 
variable  y  "'hich  is  restricted  to  four  values  which  may 


"''.■>?.'!. 
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be  taken  as   1,2,3,4,      These   fionctions    are    c?5lled    splnor 

components.      The   Kermltian   Inner   product   of   f(x,Y)    and 

4   r     

g(x,Y)  is  given  by  ;>-  /  f(x,Y)  g(x,Y)dx   and  the  sym- 

4   r 
metric  inner  product  by   I>   /  f(x,Y)  g(x,Y)dx  .   An 

Y=l  '-' 
operator  A^'*^  may  be  represented  by  an  integral  operator 

with  a  kernel  ^(x,yjx',y')  given  by  the  definition  (15), 
For  some  purposes  such  a  representation  of  t^-e  operator 
is  useful, 

Dirsc's  wave  equation  for  the  electron  in  an  elec- 
tromagnetic field  is 

(25)  i<^'{/(x  Yvt)   =  H^,Y  ^(x,Y;t) 

where 

(26)  H^'"^  =  ^  a^  \.L^   +  e  A.(x))  +  e  <l)(x)  -  mp^ 

j=l   "^    \  j      ■^      ^ 

Here  olX  and  [3^  are  Hermitian  operators  which  operate  T«'ith 
respect  to  the  y  variable  only,  Thejr  satisfy  the  follow- 
ing commutation  relations 

J   1    i   J        Ij 

(07)    \    f3^  a^  +  a*|J  p^  =  0 

Y 

where  I   is  the  identity  operator. 

The  operators  o.'C  ,  R'  can  be  expressed  as  Integral 

operators  with  kernels  ci,.(y,Y'),  S(y,y')  which  are  the 
well-known  Dirac  matrices. 


•oil  or.'! 


r.  J' 
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In  the  expression  (26)  we  have  taken  h  =  c  =  1,  The 
mass  of  the  electron  is  m  and  its  charge  is  e.  The  func- 
tions A.(x)  and  (})(x)  are  the  vector  and  scalar  potentials 
of  an  electromagnetic  field  and  are  taken  as  real, 

V'/e  shall  look  for  solutions  of  equation  (25)  which 
can  be  written  as 

(28)  vKx,Y;t)  =  e'^^*-X(x,Y;E) 
so  that  equation  (25)  leads  to 

(29)  H^'Yy(x,Y;E)  =  E  X.(x,YjE)   . 

This    equation   is   of  the   form  (18)   with   H^'"*'  (riven  by   (26), 

X  Y 

For  the  case  of  the  Dirac  electron  we  define  H  *'  as 

o 

X  Y 

and  q'  "  by 

5 
(31)      q''*^  =  e  ;^  aY  A,(x)  +  e  ^(x) 

j=l   J   ^ 

v'here  A.(x)  and  (j)(x)  are  to  vanish  for  1x1  >  r   so  that 
J  '  '  -  o 

q'^'^  =  q^ix)    =   0   for    |x|  >  r 
'Ve  have  thus  specified  H  *'  snd  q  '   of  equation  (19). 

B.   The  Incident  "Vave. 

As  in  the  general  procedure  of  Part  II  of  the  pres- 
ent note  we  shall  write 


I5'.''5T 


Y  ^  -  ■ 


r-  \ 
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"i/e  require  that  /,  4   be  a  solution  of 

where  H  »'  is  given  by  (30).   Suitable  solutions  are  the 
"spinor  plane  wave"  solutions  which  have  the  form 

(^■5^)    'Ain(x,Y;E)  =/(Y;E;7,r)  e^l^U^'^) 

where  V/  is  a  unit  vector  which  specifies  the  direction 
of  propagation,  ()'x)  is  the  inner  product  of  the  vectors 
X  and  Y/  ,    and  |k|  is  the  absolute  value  of  the  momentum 
and  is  given  by  the  relation 

(33)  \kf   =E-  -  m^   , 

The  quantities  X(Y»E,v,rO  are  functions  of  their  argu- 
ment?.  Here  T  is  a  variable  which  is  restricted  to  two 
values  which  may  be  taken  as  +1  and  -1.   The  significance 
of  ?"  is  that,  in  the  non-relativlstic  approximation,  it 
represents  the  z  component  of  the  spin. 

By  substituting  (32)  into  ths  equation  H^'V.   =  ^X. 
it  is  seen  that  the  functions   X(  Y>E,V  ■^2')  satisfy  the 
following  equation 

r  5  7 

(34)  <  E  +  |k|  ^  a*^  V   +  ^"fynl    Ari^,V;,T)    =  0      . 

For  the  purpose  of  the  present  note,  it  is  not  necessary 
to  give  an  explicit  form  for  the  functions  '/i  Y>2,V^  ,7^) ; 
these  functions  can  be  found  in  textbooks.   However,  it 
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will  be  useful  to  indicate  the  orthogonality  properties 
of  th-^  functions.   These  orthogonality  relation?  are 


(35) 


Y 


also 

(35a)    5!X(Y>E,V,'r)'X.(Y'i^,Vj,r) 

z 

+  ^'X(yv-e,v,^)  X(Y",-E,v?,'rr  =  5(y,y')  . 

0 

Prom  the  above  relations  it  can  be  shown  that  any  function 
f(Y)  of  Y  can  be  expressed  as  a  linear  combination  of  the 
functions 

"X(y-,E/^+1),  X(y-,E,v^,.1),  X(yV-E,^,-H),  and  "/j: YV-E,y^,-l) 

for  fixed  values  of  E  and  ''^  , 

C .  The  Intep-ral  Representation  of  the  Operator 
[E  ~  H^*'^]"-'-;  The  Integral  Equation  for  the 
Scattered  '"/ave. 

Through  expression  (32)  we  have  given  acceptable  func- 
tions y^  .    .      '•■'■e   wish  now  to  set  up  the  integral  equation 
for  /-o^  corresponding  to  equation  (23).   We  shall  ex- 
press  [E  -  ^-q*^    ^^   a^^  Integral  operator  such  that  ^ 
as  given  by  (23)  behaves  like  an  outgoing  spherical  wave 
for  large  values  of  |x|. 

For  this  purpose  let  us  consider  the  differential 


i  fo  r  • 
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equation 

(36)  (E  -  H^'*^)  f(x,Y)  =  k(x,Y)   , 

where  f(x,Y)  is  subject  to  the  boundary  condition  that 
for  large  values  of  |x|  f(x,Y)  is  to  behave  lil-re  an  out- 
going spherical  wave.   The  solution  of  (36)  can  be  written 

(37)  f(x,Y)  =  (E  -  K^'Y)-^  k(x,Y) 

/' 

=  >  ':    g(x,Yjx',Y')  k(x',Y')dx' 
Y'  J 

provided  k(x,Y)  dies  out  sufficiently  rapidly  for  1 prge 
value?  of  ix|.   The  second  of  equations  (37)  expresses 
the  inverse  operator  (E  -  H  *^)        ss  an  integral  operator. 
The  boundary  condition  on  f(x,Y)  will  determine  the  form 
of  the  influence  function  g(x,Y>x' ,y*  )  •   ''s  shall  then 
have  the  desired  integral  representation  for  the  inverse 
operator  (E  -  H^*"^ )""'■. 

The  influence  function  satisfies  the  following 
equation 

(38)  [E  -  H^'^]g(x,Yix»,Y')  =  6(x  -  x' )  5(y,y' )   • 

Let  us  now  consider  the  operator  [E  -  H  *M.   Then  from 

o 

t^e  commutation  rules  (87)  for  aj  ^  p'  one  has 


(39)     [IS  +  H^''^][E  -  H^'Y]  =  [E  -  H^'YjfE  +  h^,Y] 
o         o  o        o 

=  (E^  -  m^  +x;^)       =  dkl"  +  v^)   . 

In  fact  it  wfis  the  requirement  that  equation  (3  9)  hold, 
which  led  Dir«ic  to  the  commutation  laws  (27). 


{  oe,} 
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'Ve  shall  nov;  consider  a  solution  of  t'-^e  differential 
equation 

(40)  (E^  -  m^  +  V^)  3(x,x')  =  (  |kl^  +  V^)  3(x,x') 

=  5(x  -  X' )   . 

2 

(The  subscript  x  on  the  operator  'v   indicates  thst  the 

differentiations  are  to  be  carried  out  on  the  variable  x 
rather  than  x' , ) 

Any  solution  s(x,x')  of  equation  (40)  can  be  used  t® 
form  a  solution  of  equation  (38). 

Prom  (39)  and  (40)  one  has 

(41)  [E  -  H^'^ILE  +  H^'Y]s(x,x')6(y,Y')  =  5(x-x' )5( y,y' ) . 

o         o 

Hence  a  solution  g(x,Y>x',Y')  of  (38)  is 

(42)  g(x,Yjx',Y')  =  (E  +  H^'*^)  s(x,x')  5(y,Y' )   . 

This  method  of  obtaining  Influence  functions  for  the 
Dirac  operator  is  a  well-known  one.   A  solution  s  (x,x'  ) 
of  (40)  which  leads  to  a  solution  f(x,Y)  which  is  pn  out- 
going wave  is 

l|kl|x-xtl 
'«'     ^  =  -  4x|x-x-|     • 

The  influence  function  g   obtained  usins?  s   is  in  explicit 
form 


^i   ^- 
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(44) 


S^(x,Yix',Y' )  = 


JHM 


x-x' 


x-x' 


+1] 


^  ^  ai(r,Y')(x^-xJ^) 
■ft 


)  ^llkl Ix-x' 


ilkllxl  _( 

^  <   [E5(y,y') 


-  mO(Y,Y'  )  ■*-  E5(y,y'  )(     47^|x-x'|    ' 

The  functrion  g  represents  nn  outgoing  wave  because 
the  time  factor  is  e"    (see  (28)).   Thus  the  solution 
of  (36)  representing  sn   outgoing  wave  is 

(45)    f(x,Y)  =11   g„(x,Y;x',Y')  k(x',Y')dx'   . 

Y'  J       "" 

For   large   values   of    |xl,    f(x,Y)    tate  s   the   form 

5 

-  U'\    ->_  a  (y,Y')' ii 
i=l   ^         ^^ 

I  -i|k|(x',v., )  7    • 

-mn(Y,Y')]/  e        '^  k(x',Y')dx'  i-   . 

J 
Here  )..    is  a  unit  vector  with  components  V  >. -,  ^  defined  by 

(47)  X  =  |x|  V/^ 

and  is  thus  the  direction  of  observation.   It  is  thus 

seen  that  f(x,Y)  behaves  li'-B  an  outgoing  spherical  wave 

for  large  values  of  |x|,  as  desired. 

Havinp  found  the  form  of  the  inverse  operator 

[E  -  H  '']   which  satisfies  the  boundary  conditions  we 
o  -^ 

can  express   /,.    in  terms  of  X  by 
s  c 
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(48)  '/^(x,y)  =ZU^(^,Y;x',Y')q''  '"^  7{x',YMdx'   . 

The  factor  q''^  '^  7-(x',y')  in  the  inte.f^rand  means  that 
the  function  X(x,y)  which  results  from  applying  the  op- 
erator  q  *^  on  ^{x^y)    is  to  have  its  arguments  x  and  y 
replaced  by  x'  and  y'  respectively.   Equation  (49)  is 
the  desired  integral  equation  for  "/-.  . 

D.   The  Asymptotic  Equation  for  the  Scattered  ''''ave. 

For  large  values  of  |x|,  the  integral  equation  for 

"X   becomes 

^'^  »l|k|UI    ( 

(«)     \,(x,T)  =-"     4^|^|  |:^^[E5(Y,r') 

3 
-  |k|  ^  ^i(y,Y')^-l^  -  mp(Y,YM] 

f  -i|ki(x«//  )     ■  7  ■ 

J  e         ^■'-  q^  »Y  ^x',Y')dx'  r  . 

"/e  wish  nov/  to  write  (49)  in  a  form  similar  to  (24)  in 
order  that  the  variational  principle  mpy  be  used  as  de- 
scribed in  Part  II  of  this  note.   Looking  at  the  inte- 
gral in  (49)  we  see  that  the  integral  is  very  nearly  an 
inner  product  of  the  function  formed  by  the  operator 

X  Y 

q  *'  acting  on  a  spinor  plane  wave  and  the  function 
''/(x,y)»  provided  one  introduces  as  a  factor  the  function 
0((Y>E/^,f)  which  was  used  in  the  definition  of  the  spinor 

plane  wave  (see  ( 3P ) ) .   This  line  of  thought  motivates 

our  use  of  the  identity 


T.  I        '■  ^    ' 


».  «.       1 


X'JO 
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3 


(50)  [Ee(Y,YM  -  |k|  ^a.(Y,Y')7..  -in6(Y,Y')] 

1=1   ^        ^^ 

=  2  E  27(YiE,i7  ,'?:)7.(Y'>E,»r,,r)   . 
7r  ^ 

Before  showing  how  this  identity  is  used,  we  shall  verify 
it.   We  use  the  orthogonality  relation 

(33a)    ^7Sry^,ri/c)'WFTKn7U 

L 

+  S')((y-,-E,u,T)  74Y'i-E,^,r)  =  5(y,Y'  ) 
and  equation  (32),  which  we  write  as 

(51)  2  H  (y,y')  "A(Y';Ei^;r)  =  E -/( Y;E,V7,r) 

y' 

where  ^ 

i  ■   (  3 

(52)  H  (y;y')  =  -  |kl  ;^  a  (y,y')  )?,  -  m3(Y,Y')   . 

3=1   ^       '^^ 

Hence  from  (33a) 

(53)  H^(y,y')  =  >  H^(y,y")5(y",Y') 

Y^ 


=  E  I Z  /(Y;E,v^,t)-)C(Y»;E,v7,r) 


^  X(Yj-E,»7,r)  7r;^v^E7?7tT  f 


5 
The  operator  E5(y,y')  -  Ik]  ,>•  a.(Y,Y')l9.  ~  ^^(y,Y') 

J=l   ^       •  J 

may  be  written  E5(y,y' )  +  H  (y,y').   Multiplying  (33a) 
by  E  and  adding  the  resulting  equation  to  (53),  one 
obtains 

(54)   E5(y,Y')  +  H  (y,y')  =  2  E  ^  X^  Y;E,V,r  )  /(  Y' jE  ,^;,r) 
which  is  the  identity  (50). 


V  PC-. 
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Now  we  shall  show  that  we  have  /   expressed  in  the 

^sc 

form  (24)  for  large  values  of  |x|,  when  we  use  the  Iden- 
tity (50)  in  (49).   If  the  incident  wave  part  of  the  func- 
tion "X(x,y)  has  the  direction  10'  and  the  value  of  V   is 

T',  we  shall  denote  ')((x,y)  by  X^^>'''>^»^' »'^' ^  •   '^®  P^°~ 
ceed  to  define  the  "spinor  snherical  wave"  ©(x,Y>S,^,r)  ^7 

T7  ^Ijkl  |x| 

(55)  e(^,Yi3,Y^,r)  =  -  ^^-f^jY] X(YiE,>9,T)   . 

The  spinor  spherical  wave  is  analogous  to  the  spinor 
plane  wave  (3P),   Furthermore,  we  define  T(S;y^,X>"^' ,t'' ) 

by 

(56)  T(E-,»;>,r-,K}',T')  = 

Then  the  integral  equation  for  X   takes  the  form  for 
large  |x|, 

(57)  73,(x,y)  =  X^^i\x\y^^,r) 

=  Z  '5(x,Y-,^,173_,n:)  T(S-,)r;^,t^v|t,t')   . 
V 

The  scattered  wave  may  be  regarded  as  the  sum  of  two 
spinor  spherical  waves,  each  bexnp  characterized  by  a 
different  value  of  7r  .   The  function  T(S;t7  ,t,*,>?  ;cl)  may 
be  regarded  as  the  amplitude  of  the  spinor  spherical  wave 
in  the  direction  >?   and  with  "C  =    t",  when  the  incident 
spinor  plane  wave  has  as  its  direction  of  propagation  n 
and  its  value  of  T  is  T_.   We  have  therefore  obtained 
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"X   in  the  desired  form  (24),  where  T(E;i,  r;<1 '  ,r' )  is  an 
inner  product  to  which  a  variationa]  principle  can  be 
applied . 


E .   Application  of  the  Variational  Principle. 

In  order  to  show  how  the  second  variational  principle 
discussed  in  Part  I  may  be  used  to  find  the  amplitudes 
T(E^in,  ^j  V'/' ,T' )  we  need  only  show  what  quantities  are  to 
be  Identified  with  the  quantities  appearing;  in  equation 
(p).   We  carry  out  this  identification  in  accordance  with 
the  procedure  discussed  in  Part  II. 

Let  us  first  consider  T(E;^'(,T*,v^' ,T' )  as  being  the 
symmetric  inner  product  of  two  vectors.   According  to  the 
general  formalism  of  Part  II  of  this  note,  we  identify 
A^^jYj^,'" '  ,t' )  with  the  vector  y  of  equation  (9).   We  see 
that  the  vector  a'  is  to  be  identified  with 


q^''^X(Y',5:,v:,t)  e"^'^' ^^»^('^  ,  using  the  fact  that  q^»^  is 
a  Hermitian  operator. 

The  other  identifications  with  the  quantities  of 
equation  (9)  can  now  be  made  easi]y.   The  operator  K  is 
Identified  with  q^'"^  -  q^»'^[E  -  H^*'^]~^q^'^  and  the  vector 
a  with  q^*'^yjY>E,V|',1' )e^'^'  ^^*^' ^  in  accordance  with  the 
general  procedure  outlined  in  Part  II. 

In  order  to  write  the  integral  equation  corresponding 
to  tho  first  of  equations  (9)  in  a  convenient  form,  we 
shall  write  q  »'  snd  [E  -  H^»']~   as  integral  operators 
with  appropriate  kernels.   The  kernel  of  the  integral 
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(58)         q(x,Yix' ,y' )    =   e 


operator   corresponding   to  q    *'    is  'Arritten 

3 

+  e^ix)    5(x  -  x' ) 
It  is  to  be  noted  that 


a.(Y,Y'  )  A.(x)  5(x 


-  X') 


(59) 


q(x',Y';x,Y)  =  q(x,Y;x',Y')   , 


since  q  »'''  is  a  Hermitian  operator. 

The  equation  corresponding  to  the  first  of  equations 
(9)  becomes  in  this  notation 


(60) 


q(x,Y;x'  ,Y'  )  X(  Y'  ;E,T^t  ^^r  )  ^^\^\  <^^'  *^'  ^  ^^j 


q(x,Y;x",Y")  ^ 
Y 


q(x,Y;x\Y' )  -  ^ 

g^(x",Y''-,x'"  ,y"'  )q(x"*  ,y"'  ;x',YMdx"  dx'"  ? 

'X(x',Y'>E,V]',t')dx' 

The  integral  equation  (60)  is  merely  a  rewriting  of  equa- 
tion (48),  where  we  have  used  ^^.      given  by  (32).   The 
curly  brackets  contain  the  kernel  of  the  integral  operator 

To  find  the  function  A.  corresponding  to  the  vector 
y'  we  shall  v/rite  down  the  integral  equation  correspond- 
ing to  the  second  equation  (9).   The  sjTnmetric  adjoint 

X  Y 

operator  K'  *'  may  be  written  ^s    an  integral  operator 
whose  kernel  is  the  transpose  of  the  kernel  of  the  inte- 

X  Y 

gral  operator  corresponding  to  K  *'.   Accordingly,  the 
integral  equation  for  A_  is 


'  \- 
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(61)  Z  f  q(x,Yix',Y')  /(Y';S,^,Ty  e'^ '^  '  ^^' ''^^dx' 

=  2:  fU(^,Yix',Y»)  -  ^  )  q(x,Y;x",Y")  ^ 

[  gg(x",Y"Vx'",Y"')  q(x"*,Y"'ix',Y')  dx"  dx'"  (  • 

.  A(x',Y'-,5,v^,t)  dx' 
where 

(62)  gg(x,Y-,x',Y')  =  g^(x',Y'>x,Y)  . 

Just  as  the  integral  equation  (60)  for  'X(x»Y>E,'^,T) 
was  derived  from  the  differential  equation 

(63)  H^''^7v(x,Yi3,7,7:)  =  E7jx,Y-,E,v^,r) 

together  with  appropriate  boundary  conditions,  the  inte- 
gral equation  (61)  for  A (XjY'jSj^jt^   can  be  derived 
from  a  differential  equation  v/ith  certain  boundary  condi- 
tions.  We  shall  write  down  the  differential  equation  and 
boundary  conditions  and  sketch  the  derivation  of  equation 
(61)  from  these. 

The  differential  equation  satisfied  by  J\   is 

(64)  H^*'^A(x,YiE,r^,T)  =  E  A.(  x,y;E  ,'r^,r) 

so  that  j\    is  an  oigenfunction  of  H.   The  boundary  condi- 
tion on  /\  is  that  it  is  to  be  expressed  as  the  sum  of  an 

incident  wave  A.   ?3nd  a  scattered  A    such  that  A..^ 
'^  ^xn  so  -^  in 

satisfies 
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(65)      H^'Vv^^(x,Y;E,>^,r)  =  E,\^^(x,Y>E,V/,r) 

cik|x| 
nnd  A    behaves  like  the  sphorlcal  wave  « — i —  for 

large  values  of  Ix], 

We  tf5kG  PS  p  suitpbls  solution  A  .   of  (65)  the 

in 

function  'J..^   where  "/   is  given  by  (32).   The  function 

A    as  can  bo  shown  from  (64),  satisfies  the  equntion 
sc 


(66)      (E  -  hJ'Y)A^^  =  q^'^  A 


which  can  be  written  ss  the  integral  equation 


(67)      /\    =  (E  -  H^>^)"^  q^'^  A 

where  the  inverse  operator  (E  -  H  *^)         is  represented 

by  an  integral  operator  with  a  suitable  kernel.   In  a 

derivation  similar  to  that  for  g  ,  this  kernel  c^n  be 

shown  to  be  g   defined  by  (68).   Equation  (61)  is  merely 
s 

a  revolting  of  (67)  using  the  function  /\ .   described 

■*  *in 

above. 

Having  considered  the  case  where  T(E j'')  ^'rjv,' ,X' )  is 
a  symmetric  inner  product,  we  sh^ll  now  discuss  the  case 
in  which  this  amplitude  is  considered  a  Hermitian  inner 
product  of  the  two  vectors   q^»^7^(  Y*,E  ,i7.,7)  e^ '^  '  ^ '^'^^ 
and  ')(^(x,Y',E,Y  ,7:' )  .   The  identification  of  the  vectors 
a,  and  y  r:nd  of  the  operator  K  are  as  before  in  the  case 
of  the  symmetric  inner  product.   The  vector  a'  is  now 
identified  with  q^*'^ 'Yiy'yE  ^Yf^Z)  e^^^^^^'^'  'y    the  operator 
K'  which  is  the  Hermitian  conjugate  of  K  is  taken  as  an 
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integral  operator  whose  kernol  is  the  complex  conjugate 
of  the  tnnsposcd  kernel  of  the  integral  operator  which 
represents  K.   The  vector  y'  is  identified  with  the 
function  ^^(x,y,E,)),f )  which  satisfies  the  following 
integral  equation  which  corresponds  to  the  second  of 
equations  ( 9) . 

(68)    Z  fq(x,Yix',Y0  7v(Y",E,7,r)e^l^|(?'^'^dx' 
Y'  J 

=  zf  iq(x,Y>x',Y')  -  Z/q(x,Yix",Y") 
Y'j  L  r^J 

5,  I  g^(x",Y"-,x'"  ,y'"  )  q(x'"  ,y'"  Jx^yM  dx"  dx'"  i 

n(x',Y'-,E,Y,?r)dx' 


where 


(69)      g+.(x,Y;x',Y' )  =  g^U',Y';x,Y7 


r 


3 


=  ,\U^^llx-x'l-i]   -^a(YY')(x   -xM 


7   -ilkllx-x' 


-mS(Y,Y')  -E6(Y,Y')j  ^  4^1^,^,  I 


We  can  express  jD-(x,y>^  »^»f )  as  a  solution  of  a  dif- 
ferential equation  with  suitable  boundary  conditions.   We 
shall  sketch  the  method  by  whicJi  the  integral  equation 
(68)  may  be  obtained  fromi  the  differential  equation  with 
the  appropriate  boundary  conditions. 

The  function  n(x,YiE,v^,-jr)  satisfies  the  differential 
equation 

(70)        H^»YrL(x,Y>E,y,r)  =  ErL(x,Y-,E,7,r) 
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so  that  I~l  ,  like^C  >  is  an  eigenfunction  of  n    *' .      How- 
over,  the  function  n  has  different  boundary  conditions 
th^n  /  .   The  boundary  condition  on  ^-  is  that  it  should 
be  expressed  ps  the  sum  of  an  incident  pprt  r^*-.   and  a 

"concentrpting"  pnrt  ^I   .   The  function  iTL-,   satisfies 
"         con  in 

(71)         H^'^a^„(x,YiE,7^,qr)  =  Ea.^(x,Y•,EY^,^) 


and 


is  taken  to  be  /C^   rs  given  by  (3S).   The  function 

I?L    is  specified  by  the  condition  that  \}.         is  to 
con     1-         J  ^  con 

behave  like  an  inwardly  moving  spherical  wave  for  large 

values  of  Ixl.   The  equation  for  ^L    is 
'  '  con 

(72)      (E  -  H^''^)ri,on  =  ^''''—   • 
Equation  (72)  can  also  be  written 

(73)    r^,,„=  [E  -H^.r,-l,x,r-^   _ 

Using  the  boundary  condition  that  _n.    represent  an 
^  -^  con   ""^ 

inward  moving  spherical  wave  for  large  values  of  |xl, 
the  inverse  operator  [E  -  H^»']""   cnn  be  expressed  an  an 
integral  operator  with  the  function  g   defined  by  (69) 
as  kernel.   Equation  (68)  is  then  simply  another  form  of 
the  equation  (73)  usinf^  the  integral  operator  represen- 
tation of  the  inverse  operator.   It  might  be  noted  that 
g   is  the  solution  of  equation  (38)  expressed  in  the 

form  (42)  when  the  solution  s  of  (40)  is  taken  to  be 

^-i|k||x-xM 
t  ~  "   4-11 1  x-x'  I 

instead  of  s   given  by  (43), 
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